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ABSTRACT 

Explosions of massive stars are believed to be the source of a significant fraction of gamma-ray 
bursts (GRBs). If this is indeed the case, then the explosion blast wave propagates into a complex 
density structure, composed of a stellar wind bounded by two shock waves - a wind reverse shock 
and a forward shock. As the explosion blast wave reaches i?o, the radius of the wind reverse shock, it 
splits into two shock waves - a reverse and a forward shock wave. We show that the reverse shock thus 
produced is not strong, therefore full analytical treatment is required in calculating its properties. We 
calculate the dynamics of the flow and the evolution of the blast waves in all of the different stages. 
We show that the fluid Lorentz factor at r > i?o is equal to 0.725 times the blast wave Lorentz factor 
as it reaches i?Oj and is time (and r) independent as long as the blast wave reverse shock exists. 
Following the calculation of the blast wave evolution, we calculate the radiation expected in different 
energy bands. We show that about a day after the main explosion, as the blast wave reaches i?o, the 
observed afterglow flux starts to rise. It rises by a factor of about 2 in a few hours, during which the 
blast wave reverse shock exists, and then declines. We show that the power law index describing light 
curve time evolution is different at early (before the rise) and late times, and is frequency dependent. 
We present light curves in the different energy bands for this scenario. 

Subject headings: gamma rays: bursts — gamma rays: theory — plasmas — radiation mechanisms: 
non-thermal — shock waves 



1. INTRODUCTION 

In recent years, there is increasing evidence that long 
duration (igo > 2 s; Kouveliotou et al. 1993) gamma-ray 
bursts (GRBs) are associated with the deaths of massive 
stars, presumably arising from core collapse (Woosley 
1993; Levinson & Eichler 1993; MacFadyen & Woosley 
1999; MacFadyen et al. 2001; Zhang et al. 2003). 
This evidence includes the association of some GRBs 
with type Ib/c supernovae [ GRB98G425 and SN1998bw 
(Galama et al. 1998); GRB011121 and SN2001ke (Gar- 
navich et al. 2003); GRB030329 and SN2003dh (Hjorth 
et al. 2003; Stanek et al. 2003); and GRB031203 and 
SN20031w (Malesani et al. 2004)], as well as the associa- 
tion of GRBs with massive star forming regions in distant 
galaxies (Paczynski 1998; Wijers et al. 1998; Fruchter et 
al. 1999; Trentham a/. 2002). Further clues arise from 
evidence of high column densities toward GRBs, which 
associate GRBs with molecular clouds (e.g., Galama & 
Wijers 2001). 

If indeed GRBs are associated with the deaths of mas- 
sive stars, then the circumburst environment is influ- 
enced by the wind from the star. As supersonic wind 
from the star meets the interstellar medium (ISM), two 
shock waves are formed: a forward shock wave that prop- 
agates into the ISM, and a reverse shock that propagates 
into the wind (in the wind rest frame). The "wind bub- 
ble" , thus formed is composed of the unshocked wind, the 
shocked wind and the shocked ISM (Castor et al. 1975; 
Weaver et al. 1977). This complex structure of the 
circumburst environments is expected to affect the dy- 
namics of the GRB blast wave during its late (afterglow) 
evolution, hence to have an observable effect on the after- 
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glow emission (Wijers 2001; Ramirez- Ruiz et al. 2001; 
Chevalier et al. 2004; Ramirez-Ruiz et al. 2005; Eldridge 
et al. 2005). 

The main effect is expected to take place when the 
relativistic blast wave reaches the density discontinuity 
produced by the wind reverse shock. The blast wave 
then splits into two shock waves - a blast wave forward 
shock and a blast wave reverse shock (Sari & Piran 1995; 
Ramirez-Ruiz et al. 2005). The blast wave reverse shock 
thus produced, propagates into a wind that was already 
shocked by the original blast wave prior to its split. This 
wind is hot, therefore the ratio of the energy densities, 
or of the gas pressures downstream and upstream of the 
flow past the blast wave reverse shock is not much greater 
than unity, hence the blast wave reverse shock is not 
strong. Obviously, as the blast wave splits into two, the 
dynamics of the created shock waves can no longer be 
described by a self-similar motion (Blandford & McKee 
1976), which determines the evolution of the original rel- 
ativistic blast wave at earlier stages, as well as the dy- 
namics at a much later stage. 

In this paper, we analyze in detail the effect of the 
circumburst environment on the dynamics of the blast 
wave(s). Following the analysis of Castor et al. (1975), 
we determine in §2 expected circumburst conditions for 
a GRB progenitor. In §3 we use the jump conditions 
at the blast wave reverse and forward shocks to show 
that a simple analytic relation between the velocity of 
the shocked fluid prior to the blast wave split and the 
velocity of the shocked fluid after the split is obtained. 
We further find a simple analytical relation that connects 
the Lorentz factor of the reverse shock to the Lorentz 
factor of the fluid as it reaches the contact discontinuity. 
We use these results to determine in §4 the evolution of 
the fluid velocity. We calculate in §5 the resulting light 
curves in different energy bands, before summarizing and 
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discussing the implications of our findings in §6. 

2. THE CIRCUMBURST ENVIRONMENT 

The circumburst environment during the GRB explo- 
sion depends on the evolutionary stages of the progenitor, 
prior to its final (presumably, Wolf-Rayet) phase. A stan- 
dard evolutionary track for a massive galactic star is to 
start as an O star, evolve through a red supergiant (RSG) 
phase or luminous blue variable (LBV) phase, before end- 
ing as a Wolf-Rayet star (Garcia-Segura et al. 1996a, b). 
The RSG phase may be absent for low metallicity stars 
(Cliiefh et al. 2003), which are preferred as GRB progen- 
itors (Le Floc'h et al. 2003; Fynbo et al. 2003; Vreeswijk 
et al. 2004), or for rapidly rotating stars (Petrovic et 
al. 2005). 

We thus consider a massive (M > 25M0), low metal- 
licity [Z ~ 0.01) star as a GRB progenitor. During the 
Wolf-Rayet phase of the star, which lasts a duration of 
~ 10^ yr, the star loses mass at a typical mass loss rate of 
M w 10"^ A'/Qyr"'^, producing a wind characterized by 
a typical velocity ^ 1000 km s^^, presumably steady 
during most of the Wolf-Rayet phase of the star (Vink, 
de Koter & Lamers 2000; Chevaher et al. 2004; Vink & 
de Koter 2005). The evolution of the wind-driven cir- 
cumstellar shell was first derived by Castor et al. (1975) 
and Weaver et al. (1977). It was shown that during 
most of its lifetime (neglecting very short early stages), 
the system has four zones consisting, from the inside out, 
of (a) a hypersonic stellar wind with characteristic den- 
sity n„(r) = M / {AniTipr'^Vw)-, (b) a hot, almost isobaric 
region consisting of shocked stellar wind mixed with a 
small fraction of swept-up interstellar gas; (c) a thin, 
dense shell containing most of the swept-up interstellar 
gas; (d) ambient interstellar gas. 

Neglecting the width of region (c) compared to region 
(6) (see below), and assuming that most of the energy in 
region (6) is in the form of thermal energy, it was shown 
by Castor et al. (1975) that the outer termination shock 
radius is at 

rtSs,wind - (,3087r^ \ pisM ) (1) 

= 1.6x lOi^'Afl/^^.s^/Sn-J/^t^/g^cm, 

where M = 10"^M_6 Mgyr'^ Vyj = 10^w^,8cms"\ 
the ambient ISM density is pisM = rUpnisM, tiism = 
10^ no.3 cm^'' and = lO^t^.gyr is the lifetime of the 
Wolf-Rayet phase of the star. Here we have scaled the 
ambient density to a value typical of a molecular cloud, 
in which we assume the young star to be still embedded. 
The density of the swept up ISM matter in region (c) is 
approximated by its value for a strong, adiabatic shock, 
Pc ^ ApiSM- ^ Comparison of the total ISM mass swept 
up to radius Rfs,wind, Mism ^ {Att /2,)R)^.^i.^^pisM 
to the mass in region (c), « AnRlj^ ,^,-^^^ARcPc then 
leads to the conclusion that the width of region (c) is 

ARc ~ Rfs,wind/^^- 

^ Detailed models of shock propagation into ambient medium 
(e.g., Chevalier et al. 2004; Eldridge et al. 2005) suggest that 
the numerieal pre-factor depends on early stages of the stellar evo- 
lution, and may be different than 4. Nonetheless, wc show in §3 
that the explosion blast wave is not expected to reach this region 
while relativistic, and therefore the exact value of the density in 
this region has no observational consequences. 



The pressure in regions (6) and (c) is Pb = Pc = 
(2/3)uh, (assuming a monatomic gas) where Ub is the 
energy density in region (6), or 

Pb- Pc - 25 l308^ J PISM \p,sMtl ) 

= 1.4 X 10-10 m!/%^ 8^/'no4'C6^' dyncm-2. 

(2) 

The radius of the inner (reverse) shock was calculated by 
Weaver et al. (1977) (sec also Garcia-Segura & Franco 
1996), assuming that the pressure in region (6) is much 
larger than the pressure in region (a) Pb '> Pa (strong 
shock assumption), by equating the momentum flux up- 
stream and downstream of the shock. Comparison of the 
ram pressure in the upstream region (a), Pa{Ro)v'^ to the 
pressure downstream, Pb + Pbvf, leads to 

/ ■ \ 1/2 
P = R , — ( 3 Mv^ \ ' 

rcQ — ^rs^wind — I 4 ^^^p^ j 

= 1.6 X IQis M%'\^ gl/10„-3/10^2/5 

' (3) 

where PaiRo) = M /{AttRIv^), and pb = Apa, Vb = 
(strong shock assumptions) were used. The number den- 
sity of particles in region (a) is na{r) oc r~^, and at 
r = Ro is given by 

MRo) = -. ^^E2— = 3.0 X 10-2 iJ-j8M_6t;»„,8"'cm-3, 

(4) 

where Rq = lO^^iJo.is cm.* The density in region (6) 

depends on the uncertain physics of the heat conduction. 
Heat conduction could be prevented by a magnetic field, 
which is expected to be toroidal in this region (Chevalier 
et al. 2004). Under this assumption, and using the fact 
that the internal speed of sound in this region is much 
higher than the expansion velocity, the number density in 
region (b) is approximately r -independent, and is equal 
to Ub ~ 4na(i?o) (Weaver et al. 1977); (see also Dyson 
& Williams 1997). 

A schematic density profile of the bubble is shown in 
Figure 1. While being only a schematic representation, 
the density profile presented is in very good agreement 
with detailed models of stellar evolution (Chevalier et 
al. 2004; Ramirez-Ruiz et al. 2005; Eldridge et al. 2005). 

3. BLAST WAVE INTERACTION WITH THE DENSITY 
DISCONTINUITY AT Ro 

We now consider the relativistic blast wave created by 
the explosion producing the GRB. Since the rest-mass 
energy of the material in region (6), Ehm = Mti^c^ « 
2 X 10^^M_6t^,6 erg is larger than the isotropically equiva- 
lent energy released in the explosion producing the GRB, 
~ 10^^ erg, the blast wave cannot cross region (6) while 
relativistic. Therefore, all the observable effects are ex- 
pected to occur as the blast wave propagates in regions 
(a) and (6). In region (a), at radius r < Rq, the blast 
wave evolution is well approximated by the Blandford 
& McKee (1976) self-similar evolution for an explosion 
into a density gradient. A similar description holds in 
region (6) for r ^ Rq (as long as the blast wave remains 

^ Note that this equation differs from (Wijers 2001) due to an 
error in the latter. 
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Fig. 1. — Schematic density profile for the scenario of a massive 
star emitting wind. Region (a) consists of the unshocked wind, 
with a density profile n(r) oc r~^. Region (6) consists of hot, 
nearly isobaric shocked wind. Region (c) consists of the shocked 
ISM, and region (d) is the unshocked ISM. 



relativistic). We thus concentrate on the interaction of 
the blast wave with the density discontinuity at Rq. 

Consider a relativistic blast wave that propagates in 
region (a) . The matter in the downstream region of the 
shock, which we denote as region (a), is composed of the 
shocked material of region (a) . Being shocked by a rela- 
tivistic shock wave its temperature is much higher than 
nieC^, thus the pressure in this region is related to the 
energy density in the region by the relation Pa = (l/3)w5 
(as opposed to the relation used in §2 describing the flow 
in region (6), Pb — (2/3)ufc, which is valid for flow with 
temperature much smaller than m^c^, as is the case in 
region (&) prior to the blast wave propagation). Since for 
r < Rq the fluid in region (a) is downstream from the 
flow past the shock wave, it thernializes, hence isotropizes 
immediately after passing the shock. We therefore adopt 
the commonly used approximation that the fluid in re- 
gion (a) is uniform and isotropic. Under this assumption, 
its number and energy densities are na(r) ~ AT(r)na(r), 
Ua{r) ~ AT"^ {r)na{r)mp(? respectively, where T{r) is the 
Lorentz factor of the flow downstream (see figure 2, left), 
and r(r) > 1 assumed (Blandford & McKee 1976). 

As the relativistic blast wave reaches i?o, region (a) no 
longer exists, and all of its matter is swept up by the 
shock and is in region (a). At Rq the blast wave splits 
into two: a relativistic forward shock that continues to 
propagate forward into the matter in region (6), and a 
reverse shock that propagates into region (a). Thus, two 
new regions are formed: region (&) which contains matter 
from region (a) shocked by the reverse shock, and region 
(c) which contains matter from region (6) shocked by 
the forward shock. Regions ib) and (c) are separated 
by a contact discontinuity. The fluids in regions (6) and 
(c) are both at rest relative to the contact discontinuity. 
Therefore, the fluids in these two regions propagate at the 
same Lorentz factor r2 (in the observer frame), which is 
smaller than Fi - the Lorentz factor of the fluid in region 
(a) at r — Rq (figure 2, right). 

While the relativistic forward shock is strong, the re- 
verse shock thus produced cannot be strong: for a rela- 
tivistic forward shock, the energy density in region (c) 



is Uc — 4r|w;,, where ojf, ~ nbrUpC^ is the enthalpy 
in region (b)^. For a strong reverse shock, the energy 
density in region (6) is given by ui ~ (4F2 -I- 3)F2a;a. 
Here, = (4/3)4F^na(i?o)'7ipC^ is the enthalpy in re- 
gion (a), and F2 is the Lorentz factor of the fluid in 
region (6) as viewed in the rest frame of region (a), 
F2 ~ (l/2)(Fi/F2 -hFa/Fi). The term -f3 in the expres- 
sion for is added because the Lorentz factor of the 
fluid in region (6) as viewed in the rest frame of region 
(a) F2, is only mildly relativistic. Equating the energy 
densities at both sides of the contact discontinuity sepa- 
rating regions (h) and (c), using ji^ = Aua^RQ), leads to 
F2 = F2[(4/3)F2 + l]Tl. Since F2 > 1, the requirement 
F2 < Fi cannot be fulfilled. We therefore conclude that 
the reverse shock formed as the blast wave reaches Rq is 
not strong. 

The Lorentz factor r2 of the fluid in regions (6) and (c), 
the number density in region (b), n^, and the Lorentz fac- 
tor of the blast-wave reverse shock, Fj^s, are found using 
the reverse shock jump conditions and the requirement of 
pressure balance across the contact discontinuity, which 
leads to = Uc. We write the Taub adiabatic (e.g., 
Landau & Lifschitz 1959) at the reverse shock in the 
form 

_ ^j, /^ ^b - n + _ 4Fi ( V\ + i2ri \ 

nl u;,{u;a + Pb-Pa) Tj {sTi + ATj J ' 

Here, Pa = ua/3, = are the pressures in re- 

gions (a), (6) (the temperature in region (6) is higher 
than the temperature in region (a) , which is much higher 
than TOgC^, thus relativistic formulae are used), itj, = 
Uc = AT^nbrripC^ is the energy density in region (6), 
Ua ~ AT\na{RQ)mpC^ is the energy density in region 
(a), the enthalpies in regions (a), (6) are cja = wa + Pa, 
~ '^b'^-^b respectively, and we have used Ub = 4rtQ(Po) 
in the derivation of the second equality. 

By making a Lorentz transformation from the reverse 
shock rest frame to the observer frame, conservation of 
particle number flux at the reverse shock is written in 
the form 

TiaTiPi - ngr2/32 {na/ni)Ti(3i - F2/32 

PRS = 7; = — 7 1 — Tf; 7; , (0) 

J^aLi - n^F2 (n-a/n-fc)ri - F2 

where /3i.2 = (1 — ^/^i 2)^^^ is the normalized velocity of 
the fluid in regions (a),(fe) respectively, and (3fjs = (1 — 
1/F|jg)^/^ is the normalized velocity of the reverse shock 
in the observer frame. As the third equation we use the 
continuity of the energy flux at the reverse shock, which 
after Lorentz transformation to the observer frame, takes 
the form 

UaTl (1 - I3,I3rs) ilSi-pRs) = (1 - /32/3«s) (/32-/?fls) 

Equations 5, 6 and 7, which connect the thermody- 
namic properties upstream and downstream of the flow 

^ The temperature in region (6) is T;, ~ 10^?)uj_g^ K (neglecting 
radiative cooling) and may be lower if radiative cooling is consid- 
ered (Castor et al. 1975; Ramirez-Ruiz et al. 2005). Therefore 
the thermal energy in this region is much smaller then the rest- 
mass energy, and the pressure is much smaller than the relativistic 
energy density. 
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Fig. 2. — Schematic description of the different regimes during 
the blast wave evolution. Left: blast wave evolution in region 
(a) results in the creation of region (a), whose observed width 
cm is too narrow to be observed on the 
scale of the plot. Right: schematic description of the different 
plasma regimes as the blast wave propagates at Rq < r < Ri, and 
the reverse shock exists. Region (6) contains plasma from region 
(a) shocked by the reverse shock, and region (c) contains plasma 
from region (b) shocked by the forward shock. Fi > r2 > F^g. 



past the reverse shock are sufficient to calculate the un- 
linown values of the thermodynamic variables Tjis 
and rijj. The calculation is done in the following way: in- 
serting Pfis from equation 6 and using Luj^/uja = 4^2/^ ij 
equation 7 takes the form 



(f)rir2(i-/?i/32) 



1 



Equation 8 is simplified by approximating /3i_2 
l/(2r? 2), which leads to 



(8) 



1 - 



2rir2 



+ ATI 



(r? 



Tl + ATI 



Inserting the ratio na/n^ from equation 5 into equation 9, 
after some algebra we are left with a quadratic equation 
for Tl 

32ri + srlrl - isrf = o, (lo) 

with a solution 



^ 2 



\/27-l-p2 
8 -^1' 

0.725ri. 



(11) 



Inserting this result in equation 5, the number density in 
region (b) is given by 

n|^l»4zi5?.3, (12, 



l.TSria. The energy density in region (6) is 
3\/3 - 1 



4r| 



2.1, 



(13) 



which means that the energy per particle in region (6) 
is u^/n^ ~ {2.1/1.73)ua/na ~ 1.21ua/na. We therefore 
conclude that the energy per particle is increased by 
20% as the particle passes through the reverse shock, 
from region (a) to (6) . The Lorentz factor of the reverse 
shock is calculated using equation 6, 



Trs - 0.43ri. 



(14) 



4. BLAST WAVE EVOLUTION 

As long as the explosion blast wave propagates in re- 
gion (a) at r < i?o, its evolution is well described by the 
Blandford & McKee (1976) self-similar solution for an 
explosion into a density gradient n{r) oc 



„-2 



r(r;r <i?o)- 



9^; 



167rAc2 



1/2 



(15) 



where E is the (isotropically equivalent) explosion en- 
ergy, and A = M/{Attv^). At r = i?o, 



Tl EE r(r - i?o) = 20.5Eli\'/'R,',i\ (16) 



where na{Ro) = 10 ^-^ni?,, _i.5 cm ^. 

The spatial dependence of the hydrodynamic variables 
of a shocked fluid element (in region (a)), are given by the 
Blandford & McKee (1976) self similar solution. In this 
solution, as the blast wave expands to radius r, more 
than 90% of the energy and the particles are concen- 
trated in a shell of co- moving thickness Ar?°- = ^r/r(r), 
where ^ is a numerical factor, in the approximate range 
0.1 — 0.5 for the hydrodynamical quantity in question 
(number density, energy, Lorentz factor etc.). Adopting 
the approximation that the fluid in region (a) is homoge- 
neously distributed, we write the comoving width of this 
region at r = Rq as ARf {r = Rq) = Rq/SCTi, where C 
is a numerical factor of order unity, which is inserted in 
order to parametrize the discrepancy between the actual 
density and energy profiles and the homogeneity approx- 
imation used. The Lorentz factor of each fluid element 
in region (a) a,t r = Rq is therefore approximated to be 

r(r = i?o) = ri. 

At r > Rq, the flow in region (a) which was down- 
stream the flow past the blast wave at r < Rq, be- 
comes upstream the flow past the reverse shock. A 
fluid clement in region (a) therefore continues to move 
at an approximately constant (r-independent) Lorentz 
factor r = Fi, as long as the reverse shock exists ^. 
The shock jump conditions analyzed in §3 then imply 
that the Lorentz factor of a fluid element in regions 
(6) and (c) is also r-independent, given by equation 11, 

r2 = o.725ri = u.8Eli\-l^^^,,R-%\ 

The Lorentz factor of the reverse shock during its life- 
time is determined by the reverse shock jump conditions 
as well, hence is r-independent and given by equation 



14, Tns = 0.43ri = 8.8 Eli\-y\,R,^,, 
shock therefore completes its crossing through region (a) 
at distance 



-3/2 



The reverse 



i?i — Ro 



Ai?f (r = i?o) 



i?n 1 



1 



= l.G6i?o 



(17) 

from the explosion. Here, Ai??'' (r = Rq) — /S.R'^ {r — 
Rq)/Vi is the observer frame width of region (d) at r = 
i?o, and C = 1 assumed in the last equality. At r > i?i 
region (d) no longer exists as all of its content is in region 
(6), and the reverse shock disappears. 

^ The velocity of the reverse shock in the rest frame of region 
(a), Pjis = 0.69, is of course larger than the speed of sound, 
Paound = 0.57. We therefore do not expect a significant change 
in the thermodynamic properties of region (a) due to adiabatic 
expansion during the reverse shock crossing time. 
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calculations, we use the standard synchrotron emission 
model, which is in very good agreement with afterglow 
observations (e.g., Wijers et al. 1997; van Paradijs et 
al. 2000). We divide the calculation of the emitted radi- 
ation into the 3 different phases, corresponding to the 3 
phases of the blast wave evolution: (a) the earliest phase, 
which corresponds to blast wave evolution in the inner- 
most regime, r < Rq, (b) the intermediate phase, corre- 
sponding to the forward and reverse shock wave evolution 
at _Ro < < ^1, and (c) the latest phase, corresponding 
to blast wave propagation at r > 



3x 10' 



10'° 

- radius from progenitor 



3x 10" 



Fig. 3. — Schematic description of the plasma Lorentz factor as a 
function of r. For r < Rq, r(r) oc r~^/^; At Rq < r < Ri, plasma 
in region (a) continues to move at Fi = F(r = Ro), while plasma in 
regions (6), (c) moves at Lorentz factor F2 = 0.725ri. The reverse 
shock moves at Lorentz factor F^js = 0.43Fi. At r > plasma 
in regions (b), (c) moves at a self-similar motion with Lorentz factor 
Focr-3/2. 



The total mass of the matter swept up by the blast 
wave as it propagates from Rq to i?i is 



mb{r = Ri) = ^mpUb (i?? 



# -1 



0.24ma, 



(18) 



where rua — ATrnip J^^^na{r)r'^dr is the total swept up 
mass of region (a). The thermal energy of the parti- 
cles swept up as the blast wave propagates from Rq 
to Eth — T^mbC^ (assuming no radiative losses) 
is much larger than the fluids' kinetic energy at 
Ek ~ T2{ma + mb)c^ . As long as the reverse shock ex- 
ists, the energy excess is compensated by kinetic energy 
loss of particles moving through the reverse shock, from 
region (a) to region (&). At r > it!i the reverse shock 
no longer exists, and therefore the plasma decelerates. 
Since at this stage the thermal energy is already much 
larger than the kinetic energy, a self similar expansion 
follows soon after the reverse shock ceases to exist, at 
r ~ This self-similar expansion is well described 

by the Blandford & McKee (1976) solution for an ex- 
pansion into a uniform density medium. We thus con- 
clude that at r > i?i the fluid Lorentz factor is given by 
r(r > i?i) = r2(r/i?i)-3/2. 

The evolutionary stages of the fluid's Lorentz factor in 
the different regimes are summarized in figure 3. The 
initial blast wave exists at r < i?o- AX Rq < r < Ri 
fluid in region (a) moves at a constant Lorentz factor Fi, 

while fluid in regions (6) and (c) moves at a constant 
Lorentz factor At r > the fluid occupies regions 
(6) and (c) only, and moves at a self similar motion. The 
blast wave becomes non-relativistic (r2 — 1 ~ 1) at radius 
Rnr ^ 3.3i?i ~ 3.5 X lO^s i?o.i8 cm. 



5. AFTERGLOW RE-BRIGHTENING 



The complex dynamic of the blast wave evolution has 
observational consequences. In this section, we calcu- 
late the expected light curve in this scenario. In our 



5.1. Emission during the early phase, r < Rq 

The Lorentz factor of the shocked plasma (in region 
(a)) at r < Rq is given by equation 15. The character- 
istic time at which radiation emitted by shocked plasma 
at radius r is observed by a distant observer was calcu- 
lated by Waxman (1997) for the case of an explosion 



into a uniform medium, t 



ob. 



r/AT c. Repeating the 



Waxman (1997) calculation for the case of an explosion 
into a density gradient n{r) cx r~^, we show in §A that in 
this case this relation is slightly modified, t°^' ~ r/2T^c. 
Radiation observed at time t°^- is therefore emitted as 
the blast wave approaches radius r = (9i?i°^ /87rylc)^/^. 
Radiation emitted at r = Rq is observed at time 
SttAc 



lOb. 

To — 



9E 



-R'q ~ 4.08 X 10*£;53Vi?o,-1.5i?o,18 s, (19) 



or about 0.47 days after the explosion. 

Denoting by and cb the fractions of post-shock ther- 
mal energy density, Uint = 4T'^ {r)na{r)mpC^ that are car- 
ried by the electrons and the magnetic field respectively, 
the magnetic field (in the fluid frame) is given by 

1/2 



= 7.7 X 



9E \ 

ilT I 

10-2 



( SttAc 



3/2 



P.-1/4 1/2 Ob. 



-3/4 3/4 



1.5^ 



?3/2 p 
•■0.18 ^' 

(20) 

where t"^- = lt°l'g^y day. The characteristic Lorentz fac- 
tor of the shock wave accelerated electrons is ^char — 
^e{mp/me)T{r){{p — 2) / {p — 1)), where p is the power 
law index of the accelerated electron energy distribution, 
rig (7) oc for 7 > "fchar- The resulting synchrotron 
emission peaks at 

iy^ {t°''- < ig'' ) = 5.6 X 10^2(1 -I- z)-i 



^^53 ^e. 



.1/2 



fob. 

2''day 



-3/2 - 



(21) 



i^B 'i''-^--' Hz, 

characteristic values eg = 
2 are taken (e.g., Wijers & 



where z is the redshift 
lO-ifie,-!, en = lO'^eB. 

Galama 1999), and power law index p — 2.5 assumed 
This frequency is below the break frequency of the spec- 
trum, corresponding to emission from electrons for which 
the synchrotron cooling time is equal to the dynamical 
time, tdyn r/8CcT{r), 

L'f {t°''- <t°'' ) = 7.2 X 10"(l-|-z)-i 

-3/2>2,o6. l/2„-2 D-4 



X-^53 ^B, 



2^0 ''day 



1.5-^0,18 

(22) 



The synchrotron self absorption frequency is 



,ob. (j-oh. 



{t""- < t^"-) ^ 7.8 X 10**(1 + z)- 



^''day 



1/5 ,-3/5 



-fto,- 



S,-2'50 



pl2/5 TT 

1.5-"-0,18 



(23) 
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The maximum Lorentz factor of the accelerated electrons 
is given by equating the synchrotron loss time to the 
acceleration time, tacc — E/{cqB). Synchrotron emission 
from these electrons peaks at 

i^^Ut"''- <tf ) = 9.6 X 1023(1 +^z)-i 

^-^53 ^day '^flo, -1-5-^0,18 

(24) 

Finally, the observed specific flux at z/°f- is given by 
the number of radiating electrons Ntot{t) = '^'nAr/mp 
times the observed power (per unit frequency, at v^') of 
a single electron T{r)\/^q^ B / rUeC^ , 

(25) 

where (11 = 10 dL,28Cia is the luminosity distance to 
the GRB. 

5.2. Emission in the intermediate phase, Rq < r < Ri 

At to'' the blast wave reaches Rq, and splits into the 
relativistic forward and reverse shock waves. As we 
showed in §4, the Lorentz factor of the plasma is r- 

independent in this regime. We show in appendix A, 
that in this case the time delay suffered by photons emit- 
ted at radius r compared to photons emitted at r = is 
At"''- « r/r^{r)c. We thus find that this phase lasts a 
duration 

At"''- = tf- - tl''- = 



Tic 



= 8.58 X 103 £;5-3iCo"'nij„,_i.5i?li8S, 

(26) 

or about 0.10 days. 

For Rq < r < Ri there are 3 distinctive emitting re- 
gions: particles in region (a) that were the only source 
of emission at r < Rq, continue to emit; particles in 
region (6) which crossed the reverse shock and were re- 
heated by it; particles in region (c) that are heated by 
the forward shock. All these regions are characterized 
by different thermodynamic quantities, and therefore the 
emission pattern in each region is unique. We thus cal- 
culate separately emission from the different regions. 

5.2.1. Emission from particles in region (a) 

Region (a), which was downstream of the flow past 
the blast wave at r < Rq, becomes upstream of the flow 
past the reverse shock at r > i?o. For r > Rq, the fluid 
can therefore only leave this region, high energy particles 
are no longer injected into the region, and no informa- 
tion enters it. The thermodynamic properties of the flow 
(e.g., energy density, magnetic fleld etc.,) "freeze out" at 
their values at r = Rq. The values of the peak frequency 
and of the spectral break frequency are therefore equal 
to their values at t"''- = t^''-, 



oh. 



< tf-) 



= 1.8 X W^{l + z)- 

X F2 ,2 ,1/2 



-3/2 



(27) 



'°''-(a;tg''- < t 



oh. 



< tf-) 



4.9 X 101^1 + ^) 



xe 



3/2 >2^-3/2 
2*50 



^0, 



1.5-^0,18 



(28) 



The comoving width of region (5) decreases linearly from 
its width at to, Ai?e° (t = to) = critg''74C to at tf-. 
Therefore, the optical depth decreases linearly with time, 
and the synchrotron self absorption frequency is given by 

7 



'""Lia-^tf' < < if ■) = 1-2 X 10' 



xCo 



1/5 1 1/5 
'53 ^6,-1^3,-2 
3/5 3/5 



3/5 

IT, (1+^) 



"iJo,-i.5Hz- 



(29) 

Since energetic electrons arc not injected into this 
region at t°''- > tg''-, the highest energy electrons, 
having Lorentz factor 7max(i°''' = tf-) = 3.2 x 

10* Sgg^/^Cs^^t-Roas = ^o*" ' cool by synchrotron 

radiation. Synchrotron emission from these electrons 

peaks at 



<''ax(a;C < if) 

= min [u^.^it"''- = tf-) 



l.._tol..12 (l + ^) 



'53 ^ 
1/2 



vF-2,-3/2a2 
-^-^53 ^S,-2'»0 



^ '^flo, -1.5-^0,18 Hz 



where iy^^^{t' 



z)-^E, 



l/2„-l/2 
53 ^Ra 



j-ob. 



1.5-^0 18^ Hz is the frequency of the 



(30) 



maxv ''0 J 



1.1 X 10^4(1 



synchrotron emitted photons from the most energetic 
electrons at f''- ^ tf-. At f''- = tf-. 



^S?ax(«;i°'- = tf -) = 6.0 X 1019(1 + z)-i 

^^B,-2^0 '''_Ro, -1.5^0, 18 1^^- 

(31) 

The reverse shock crosses this region at constant velocity, 
therefore the number of radiating electrons in this region 
decreases linearly with time. As a result, the total flux 

ob. 
' 



j-Ob 
j.ob. 



at i'!^- decreases linearly from its value at t"''- = 



Fm^^{a-X''- = tf-) = IS.BE^se'^ 
at t"''- = tf-. 



2"L,28"'Jio,-1.5 



mJy to 



5.2.2. Emission from particles in region (b) 

Assuming that similar fractions eg = 10^^ of the post 
shock thermal energy are converted to magnetic fleld at 
both the reverse and the forward shock waves, the mag- 
netic fleld produced by the reverse shock wave 



Bg(prod.) = {STTeBnit"''- > tf-))'/^ 

= (STres X 2.lMa(t°^- > 0)'^' 



(32) 



(see eq. 13), is smaller than the magnetic fleld advected 
with the flow from region (a). 



B-,{adv.) = ^Bait""- > tf-) 

= V^Bait"''- > tf-) 

,,1/2 1/2 p-3/2 , 
^53 tB,-2-"'0,18 



(33) 



by a factor ^3/2.1 « 1.2. Here, the energy density and 
the magnetic field in region (a) assume their value at tf- 

The photons have to cross several layers of plasma (regions (6) 
and {(■')) before reaching the observer. The observed self-absorption 
frequency therefore depends on the plasma parameters at the dif- 
ferent regions. 
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for t""- > t^"-. Wc therefore assume that the magnetic 
field in this region is equal to the advected magnetic field, 
given by equation 33. 

The re-heating of the plasma by the reverse shock 
implies that the characteristic Lorentz factor of 
the electrons in this region is larger by a factor 
{ui/ni)/{ua/na) = 2.1/1.73 = 1.21 than the characteris- 
tic Lorentz factor of electrons in region (a). The charac- 
teristic synchrotron emission thus peaks at 

< t°^- < tf ) = 3.2 X 10"(1 + z)'^ 

^El,el_,e]^% (34) 

^ "ilo,-l.5-'^0,l8 
The reverse shock re-randomizes the energy, thus the 
spectral break frequency is independent on the break fre- 
quency in region (a), and is given by 



8.6x10 



xE, 



-2-3/2^2 



^53 "^5,-2*; 
XTljj 1 K^in 18 Hz. 



(35) 



''Ro,-1.5-^a 



At t"^- = tl^-, the comoving width of this region is 0, 
hence the break frequency vf'-{h]t'^- t°^') — > oo. At 
^ob. _ ^ob. ^ ^]^-g frequency is equal to 



yob-(p-t°^- =tf-) = 1.2 X 10i9(l -hz)-i 



xe. 



3/2^4 -3/2 



flo,-1.5-"-0,18 



Rn 1 o Hz. 



(36) 



The self absorption frequency is 

ob.lh.+ob. ^ +ob. ^ +ob.\ _ ^.1 X 10^ (t°'' 

x(l + z 



t'r,(6;i8 < t < *i ) 



j.ob. 
^0 



.3/5 



p.4/5 1 1/5 
^53 ^6,-1^3 -2 



(37) 

while the highest frequency of the synchrotron emitted 
photons is expected at 



5.5 X 10^3(1 -Hz)-i 



1/2 -1/2 



53 



'^iJo!-l. 5-^0,18 I^^- 



(38) 

The observed flux at v^' emitted from this region is 
expected to grow linearly from at t^^- to 

= V3 X 0.725 X Fm^^a; t°^- = t^^-) 

= 23.2 £'53eB,-2f^L,28'^i/o%i.5 "^^^ 

(39) 

at t°^- = tf -. 

5.2.3. Emission from particles in region (c) 
While the energy density in region (c) is equal to 
the energy density in region (6), there is no advected 
magnetic field term. Hence, the magnetic field in 
this region is assumed to be produced by the forward 
shock, and is equal to J3g = Bi{adv.)/1.2 = 1.9 x 

10-i£;5^^e^(i2^o;if G (see eq. 33). Electrons are ac- 
celerated by the forward shock to characteristic Lorentz 
factor -fcharic) ^ ee{mp/me)T2{{p - 1) / {p - 1)) , produc- 
ing synchrotron radiation at characteristic frequency 

v^-{c]tf- <t°^- <tf>-) = 1.0 X 10"(1 -^2)-i 



X-E53ee,-ie 



1/2 



l'=S,-2 
X^iJo, -1.5-^0,18 



(40) 



The characteristic spectral break frequency is v°''-{c) = 
vf-(lb) X {B-JB,f, or 



vf-{t,tf- <t"^- <tf-) = 



which is equal to 



1.5x10^ 



z 



(t<"'--tS''-)^' 

V ^-2^-3/2^2^1/2 p6 TT^ 
-^-^^53 ^3,-2*50 '* fio,-l-5-'^0, 18 

(41) 



yob-{c-t°^- =tf-) = 2.0 X 10^^(1 + z)-i 



3/2 j-y ^ r I 

-^^B,-2S,o"-_Ro, -1.5^0,18 



(42) 



at t°^- = tf- . The synchrotron self absorption frequency 
is given by 



2.5 X 10^ {t"^ 

x(l + 



3/5 p-12/5 



-l"^B, 



xCo" ' ^1 



•0,18 



Hz, 



(43) 

and the highest frequency of the synchrotron emitted 
photons is independent of the magnetic field strength, 
and therefore equals to the highest frequency of the syn- 
chrotron emitted photons in region (6), v'^^JyC\t'Q - < 
fob. < ^06.-) ^ i^^^^{b;t^''- < t°^- < tf-) (see eq. 38). 

The number of particles swept by the forward 
shock from region (5) into this region grows linearly 

with time, and is equal to A'2 = 47r J^^ r'^ribdr = 

(4/3)iVi[(i?i/i?o)^ - 1] « 0.247ViCo"^ ai r = Ri (see eq. 
18). Here, A^i = AirARo/mp is the total number of par- 
ticles in region (a), that were swept by the blast wave 
at r = Rq. The observed flux at u^- emitted from this 
region therefore grows linearly with time, from at t^"- 
to 

2<^Z,28'^flo!-1.5mJy' 

(44) 



^. j-ob. 



F (r- f 
at f''- = tf 



tf ) = 4.8 E53e]^\dvl,n 



5.3. Emission in the latest phase, r > R\ 

As explained in §4, the blast wave evolution at r > i?i 
is approximated by the self-similar evolution, T{r > 
i?i) = r2(r/i?i)-3/2. At r = the shocked shefl con- 
tains matter shocked by the reverse shock in region (b), 
as well as matter shocked by the forward shock in region 
(c). While for r > Ri hot matter continues to enter re- 
gion (c) through the forward shock wave, matter in region 

(b) cools adiabatically. Therefore, the energy densities of 
the flow at the two sides of the contact discontinuity that 
separates regions (6) and (c) are different for r > Ri. For 
an explosion into a constant density medium, the energy 
density in region (c) evolves as Uc oc r2 (r)n(r) oc r~^ , 
while the energy density in region (b) is ?ig oc r"**. It 
thus follows that the contact discontinuity can not sepa- 
rate these regions any more. Matter in these two regions 
starts to mix at r > Ri, and region (b) eventually disap- 
pears. 

We thus treat the matter in regions (6) and (c) as 
having similar thermodynamic properties at r > Ri. 
This mater is concentrated in a shell, who's observed 
width at r = Ri is Ai?(r = i?i) = cAt°'' (/32 - Prs) = 
0.20 X i?i/(4r2). During the self similar motion that 
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follows, the width of the shell is AR cx We there- 
fore expect that the shells' width increases to its termi- 
nal value AR{r) = r/(4r^(r)) by sweeping matter into 
region (c), shortly after tf''. We therefore neglect emis- 
sion from region (b) compared to emission from region 
(c) at r > iii, as well as the numerical factor 0.20, and 
calculate the emission in the latest evolutionary phase 
in accordance to the self-similar solution, which is the 
asymptotic solution for r » 

The time delay suffered by photons emitted in this re- 
gion compared to photons emitted at r = is therefore 
At°^- pa r/4r^(r)c, hence the time delay compared to 
photons emitted at Ri is At°''- - iii/(4r|c). These pho- 
tons are therefore seen at time 



t 



ob. ^ob. 



Ri 



+ 



(45) 



^ 4Tlc ' 4r2(r)c 
from the explosion. This equation is written in the form 



fob. = 



4r2c' 



(46) 



where t°'>- = f''- - tf- + Ri/{ATlc). 
Using this scaling, the magnetic field is given by 



,1/2 



(47) 

where t"^- = It^^^ day. Synchrotron radiation therefore 
peaks at 



v^ {t°^- > tf-) = 3.5 X 10i2(l + z)-^ 

j,l/2 2 1/2 .ol 
^53 ^e,-l^B,-2''daj/ 

The break frequency in the spectrum is at 
j,ob.(j.ob. > ^ob.) = 3 2 X 1016(1 + ^;)-i 

xR-V2,-3/2.2.o6. TT7 

(49) 

and the self absorption frequency is 
^'ssa(^°''' > tf ) = 3-5 X 10^(1 + z)-i 

(50) 

The highest frequency of the synchrotron emitted pho- 
tons is expected at 

<ax(i°''- > tf ) = 6.5 X 1023(1 + z)-i 

The total number of radiating electrons at radius r » 
Ri is approximated by N{r) = (4/3)iVi [(r/i?o)^ - 1] « 
(4/3)A^i(r/i?o)'^, thus the asymptotic value of the flux 
emitted at v^' is given by 

i^max(i°''- > tf ) = 38.7 E,se]i%dll,n'£^_, ,mJy 

(52) 

(a somewhat lower value is expected at r > with no 
simple dependence on the values of the unknown param- 
eters). 



5.4. Predicted afterglow light curves 

Light curves are calculated using standard formulae for 
the flux at a given frequency (e.g., Sari et al. 1998), 



F„ 



.{vc/ym)-^^-^^'''{y/yo 



-p/2 



^ssa < < I'm, 
Vm<l' < fc, 
Vc<y < fmax- 

(53) 

The calculation is done separately for each of the 3 differ- 
ent phases. In the intermediate phase t°^- < t°^- < tf'', 
emission from the 3 different emitting regions is calcu- 
lated separately and summed up. 

The resulting light curves for radio (10^ Hz), mid-IR 
(10^3 Hz), optical (IQi'"^ Hz) and X-ray (lO^* Hz) frequen- 
cies axe presented in Figure 4. In producing this flgure, 
characteristic values of i?o,i8 = 1-6, f^i?o,-i-5 = 0-4 were 
assumed, which result in t^^- w 0.6 day and At"*"' = 
0.13 day (see eqs. 3,4, 19, 26). The radio band is at the 
characteristic frequency of the synchrotron self absorp- 
tion frequency at all times, the mid-IR frequency is close 
to the peak emission, the optical band is typically higher 
than • and lower than , while the X-ray frequency 
is above v"^- for the characteristic parameters chosen. 

Two important features of this scenario are seen in 
this figure. The first, which is common to all frequen- 
cies, is the increase of the flux by a factor > 2 from t^^' 
to tf'- caused by the simultaneous radiation from the 
3 different regions during this period. The second fea- 
ture is the flattening of the power law index a of the 
flux time dependence (Fi, cx t^", where a = a{iy) is 
frequency dependent) from its value at t°^- < t^^- to a 
value smaUer by 1/2 at late times, t°^- > tf-. This flat- 
tening is caused by the change of the ambient density 
profile, from n{r) cx r"^ at t°^- < f^- to constant density 
profile, n{r) cx r^ at late times of the blast wave evo- 
lution. The power low index therefore changes accord- 
ingly: for V < Vssa, Fv oc t^ if n(r) cx r~'^, and F^, oc t^/"^ 
if n{r) cx r*^; for Vgsa < v < i^m, F^ cx if n{r) cx r~'^ , 
and oc t~^/'^ if n(r) oc r°; and for f ^ < < Vci 
F^ oc if n(r) oc r'^, and F^ oc i-(3p-3)/4 if 

n{r) cx r'^ . At high frequencies, the light curve time de- 
pendence is not changed by the change of the density 
profile, F^ oc t-(3p-2)/4 foj. > 

6. SUMMARY AND DISCUSSION 

In this work, we considered the scenario of a massive 
star as a GRB progenitor. This scenario results in a com- 
plex structure of the circumburst medium, which is com- 
posed of 4 different regions: unshocked wind, shocked 
wind, shocked ISM and unshocked ISM. We showed that 
the main effect takes place when the blast wave reaches 
the wind-reverse shock which separates the unshocked 
wind and the shocked wind regions. The blast wave 
then splits, and a blast wave reverse shock is formed. 
We showed in §3 that the blast wave reverse shock, 
which separates two hot regions is not strong. By solv- 
ing the equations describing the shock jump conditions, 
we showed that as long as the reverse shock exists, the 
shocked plasma moves at constant velocity, and simple 
analytic relation between the Lorentz factors of the flow 
at the different regimes exists (eq. 11). We then cal- 
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Fig. 4. — Schematic description of the hght curves at different 
frequencies, Radio (10^ Hz), mid-IR (lO" Hz), Optic (lO^^ Hz) 
and X-ray (10^* Hz). Power law index p = 2.5 of the accelerated 
electrons assumed, and luminosity distance dj^ = 10^* cm assumed 
in producing this plot. 



culated the blast wave evolution, and the resulting light 
curves at different frequencies (§5, fig. 4). 

The resulting light curves are significantly differ- 
ent than "standard" afterglow light curves calculations, 
which assume explosion into uniform medium, or into a 
density gradient. The resulting light curve in the sce- 
nario presented here is similar to the light curve of an 
explosion into a density gradient at early times before 
the flux rise, and to light curves obtained for an explosion 
into a uniform medium at late times. It has an impor- 
tant additional feature: during a short transition, lasting 
At°'^- /f^- ~ 1/5 of the transition time, which is expected 
to take place at ~ 1 day after the GRB explosion, the 
flux rises by a factor of > 2 at all frequencies. Such a rise 
in the flux is a prediction of this model. This model may 
therefore provide a natural explanation to the rise in the 
various optical bands fluxes of GRB030329 by a factor of 
4 observed after '--^ 1.4 days (Lipkin et al. 2004). 



The occurrence time of the transition f^' , is the time 
at which radiation that was emitted as the blast wave 
reached the radius of the wind reverse shock i?o is ob- 
served. This time is weakly constrained, being strongly 
dependent on the unknown values of i?o and of the wind 
density at this radius, n^jg . Therefore characteristic tran- 
sition occurrence times lasting between few hours to sev- 
eral days may be expected. The transition time Ai°^- on 
the other hand, have the same parametric dependence as 
f^' . We therefore expect longer rise of the flux in bursts 
for which the transition occurs at late times, compared 
to bursts for which the transition occurs at earlier times. 

The purely constrained value of the transition time 
implies that observed afterglow light curves from differ- 
ent bursts may be interpreted as resulting from explosion 
into uniform medium, or into density gradient, if only 
part of the data is available. Thus, we find that both 
interpretations (e.g.. Chevalier & Li 2000) may be con- 
sistent with the data, with preferably parameters that fit 
explosion into uniform medium at very late times. 

We used here a highly simplified model to describe the 
density profile of the ambient matter. We did not con- 
sider radiative cooling, which is highly uncertain and, if 
exists, lowers the temperature of the shocked wind gas in 
region (5), hence rises its density. We also did not con- 
sider inhomogeneities inside the different regions, which 
may also have observational consequences. The calcu- 
lation presented in §3 for the interaction of a relativis- 
tic blast wave with a density discontinuity caused by a 
strong shock is, however, general, and can be used in the 
context of supernovae calculations as well. We thus con- 
clude that the main findings of this work - the afterglow 
re-brightening by a factor of > 2, and the different light 
curves slopes at early and late times, remain valid also 
in a more complex and realistic ambient medium profile. 

AP wishes to thank James Miller- Jones and Peter 
Meszaros for useful discussions. This research was sup- 
ported by NWO grant 639.043.302 to RW and by the EU 
under RTN grant HPRN-CT-2002-00294. 



APPENDIX 

RELATION BETWEEN THE OBSERVED TIME AND EMISSION RADIUS IN RELATIVISTIC 

FIREBALL 

Waxman (1997) calculated the relation between the observed time delay At"''- of photons emitted by synchrotron 
radiation from a shell expanding relativistically and self-similarly in a uniform medium as the shell's front reaches radius 
r, with respect to photons emitted from the center of the explosion at r = 0, and found the relation At"**- ~ r'/(4r^(r)c). 
Here, we repeat the calculation for the case of a relativistic expansion in a density gradient, n(r) cx r^^. In this case, 
assuming adiabatic expansion, the Lorentz factor of the flow scales with the radius as Tlr) oc r"^^"^)/^. The magnetic 
field scales as B cx r(r)n(r')^/^ oc r~^/^, and the characteristic Lorentz factor of the electrons, ^char T{r). The peak 
frequency of the synchrotron emitted photons (in the comoving frame) therefore scales as Vm oc 'y'^har^ ^ 7-P-9/2 

The number of photons emitted by a single electron at a unit comoving time Nj is proportional to the magnetic field 
B. The total number of electrons swept by the shell (and, presumably, emit) at radius r is Ne{r) oc r^~^. Since the 
comoving time during which the plasma propagates a distance dr is df^"' ~ dr/T{r)c, the number of photons emitted 
as the plasma expands a distance dr scales as dNj/dr oc NeNj{dt'^°' /dr) oc r'^-^p/^. The fraction of photons emitted 
at frequency i^m that are observed at frequency is given by eq. 5 of Waxman (1997), df /du'/^- ~ {2Tvm)~'^ , where 
F 3> 1 assumed. Thus, the number of photons produced by the fireball shell at radius r with frequency in the range 
is 



dN^ 



oc 



F(r) 



-^vZ} oc r3(3-p). 



(Al) 



dv^- dr dr 

The delay in observed time of a photon emitted from the edge of the shell at frequency Vm and observed at frequency 
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Fig. A5. — Normalized observed flux at observed frequency Um as a function of time, for various density profile considered. The solid 
line shows the result for a density profile n(r) oc as expected in the early stages of the blast wave evolution, at r < Rq. The dash 
line presents the results for an expansion in homogeneous medium, and is similar to the results of (Waxman 1997). The dash-dotted line 
represents a scenario of r-independent Lorentz factor and magnetic field, as expected for Rq < r < Ri . 



,ob. 



with respect to photons emitted on the hne of sight is given by (for F ^ 1) 



A«.K;,.)^j5j(r^-i)«r«-(r— ""-1/2) 



(A2) 



An additional delay exists between photons emitted at radius r on the line of sight, compared to photons emitted from 
the center of the explosion at r = 0, 



r 



dv T 



Assuming that photons are emitted uniformly from a shell of finite thickness ^r/F^c, the arrival time of photons 
emitted at radius r and observed with frequency v"^' are uniformly distributed over the range t{v'^\ r) = /^tg{v°^' , r) + 
Mr{r), and t{v^-,r) + r^, where = ^r/F^c. 

Figure A5 shows the results of a numerical calculation of the flux as a function of time (compare with fig. 1 of 
Waxman 1997). The width of the shell considered is ^ = 1/8. The three plots represent the three cases considered 
in this work. Explosion into density gradient, n{r) cx r^^ is presented by the solid curve, and explosion into uniform 
medium is presented by the dashed curve. The Dash-dotted curve represents the special scenario that exists in the 
intermediate phase of the blast wave evolution at ii!o < r < Ri^ during which the Lorentz factor and the magnetic 
field are kept constant as long as the reverse shock exists. 

We thus conclude that while for an expansion into uniform medium, Waxman (1997) relation, t°^- ~ r/AT^c holds, for 
an explosion into density gradient, n(r) cx this relation is modified, and is t°^' ~ r/2F^c, while in the intermediate 
phase of the blast wave evolution the correct relation is t°^^ ~ r/F^c. 



REFERENCES 



Blandford, R.D., & McKee, C.F. 1976, Phys. Fluids, 19, 1130 



Castor, J., McCray R., & Weaver, R. 1975, ApJ, 200, L107 



Signature of wind RS on GRB afterglow 



11 



Chevalier, R.A., Li, Z.-Y., & Fransson, C. 2004, ApJ, 606, 369 
Chevalier, R.A., & Li, Z.-Y. 2000, ApJ, 536, 195 
Chieffi, A., et al. , 2003, MNRAS, 345, 111 

Dyson, J.E., & Williams, D.A. 1997, The Physics of Interstellar 

Medium (lOP: Bristol) 
Eldridge, J.J., Genet, F., Daigne, F., & Mochkovitch, R. 2005, 

preprint (astro-ph/0509749) 
Pruchter, A.S., et al. 1999, ApJ, 519, L13 
Fynbo, J.P.U., et al. 2003, A&A, 406, L63 
Galama, T., et al. 1998, Nature, 395, 670 
Galama, T., & Wijers, R.A.M.J. 2001, ApJ, 549, L209 
Garci'a-Segura, G., & Franco, J. 1996, ApJ, 469, 179 
Garci'a-Segura, G., Langer, N., & Mac Low, M.-M. 1996, A&A, 

316, 133 

Garci'a-Segura, G., Mac Low, M.-M., & Langer, N. 1996, A&A, 
305, 229 

Garnavich, P.M., et al. 2003, ApJ, 582, 924 
Hjorth, J., et al. 2003, Nature, 423, 847 
Kouveliotou, C., et al. 1993, ApJ, 413, LlOl 

Landau, L. D., & Lifschitz, E.M. 1959, Fluid Mechancs (Addison- 

Wesley, Reading, Mass.) 
Le Floc'h, E., et al. 2003, A&A, 400, 499 
Levinson, A., & Eichler, D. 1993, ApJ, 418, 386 
Lipkin, Y.M. et al. , 2004, ApJ, 606, 381 
MacFadyen, A.I., & Woosley, S.E. 1999, ApJ, 524, 262 
MacFadyen, A.L, Woosley, S.E., & Heger, A. 2001, ApJ, 550, 410 
Malesani, D., et al. 2004, ApJ, 609, L5 
Paczynski, B. 1998, ApJ, 494. L45 

Petrovic, J., Langer, N., Yoon, S.-C., & Hager, A. 2005, A&A, 425, 
217 



Rainirez-Ruiz, E., Dray, L.M., Madau, P., & Tout, C.A. 2001, 

MNRAS, 327, 829 
Ramirez-Ruiz, E., Garcia-Segura, G., Salmonson, J.D., & Perez- 

Reiidon, B. 2005, ApJ, 631, 435 
Sari, R., & Piran, T. 1995, ApJ, 455, L143 
Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 
Stanek, K. Z., et al. 2003, ApJ, 591, L17 

Trentham, N., Ramirez-Ruiz, E., & Blain, A.W. 2002, MNRAS, 
334, 983 

van Paradijs, J., KouveUotou, C, & Wijers, R.A.M.J. 2000, 

ARA&A, 38, 379 
Vink, J.S., & de Koter, A. 2005, preprint (astro-ph/0507352) 
Vink, J.S., de Koter, A., & Lamers, H.J.G.L.M. 2000, A&A, 362, 

295 

Vreeswijk, P.M., et al. 2004, A&A, 419, 927 
Waxman, E. 1997, ApJ, 491, L19 

Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, M. 1977, 

ApJ, 218, 377 
Wijers, R.A.M.J., et al. 1998, MNRAS, 294, L13 
Wijers, R.A.M.J. 2001, in Gamma Ray Bursts in the Afterglow 

Era, ed. E. Costa, Frontera. F. and Hjorth J. (Springer: Berlin) 

306 

Wijers, R.A.M.J., & Galama, T.J. 1999, ApJ, 523, 177 

Wijers, R.A.M.J., Rees, M.J., & Meszaros 1997, MNRAS, 288, L51 

Woosley, S.E. 1993, ApJ, 405, 273 

Zhang, W., Woosley, S.E., & MacFadyen, A.I. 2003, ApJ, 586, 356 



